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For Poisson input to a single server, it is shown that the stationary delay 
distribution function for last-come first-served service is equal to the dis- 
tribution function for the busy -period (the interval of time during which the 
server is continuously busy) only for exponential distribution of service 
time. A similar argument shows (hat the identity persists for a group of fully 
accessible servers, each with exponential distribution of service times — a 
result which has been a curiosity for some time. Finally, the delay distribu- 
tion for last-come first-served service and a single server with constant service 
time is derived. 

I. INTRODUCTION 

Delays for last-come first-served service were first considered by 
Vaulot for a system with Poisson input to a group of fully accessible 
servers, each with exponential service time. For this order of service, 
an arrival which is not served immediately, following the biblical edict, 
goes to the head of the waiting line. Its consideration has a natural 
theoretical interest, because, as the opposite of first -come first-served, it 
seems to be a bound for the gamut of possible service assignments, or at 
least of those with simple structure. Indeed Vaulot 1 has used it to find the 
envelope of delay functions for all service assignments. 

Very briefly, Vaulot 's formulation is as follows. Let v„(t) be the prob- 
ability that a waiting demand for service, which at a given moment has 
just become n + 1 in line, waits at least t [v n (t) is the complement of a 
distribution function; v (t) is the complement of the conditional delay 
distribution function]. Then, if a is the arrival rate, and b the service 
rate for each of the c servers, the set of differential recurrence relations 
for the v n (t) is 

Vn'(t) = bev^it) - (a + bc)v n (t) + av ll+1 (t), n = 0, 1, • • ■ . 

785 



786 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1901 

Vaulot's solution of these equations will be given later. Here it is suffi- 
cient to notice that the same relations had already appeared in the 
formulation of the busy period, for the same traffic system, by Palm. 
The correspondent to *>„(/) is /„(/), the probability that a busy period 
(all servers busy), which at a given moment has n waiting customers, 
continues as a busy period at least t. Since the boundary conditions also 
agree, v n (0 = f»(0, and in particular v (t) = /<,(/); that is, the condi- 
tional delay distribution function for last-come first-served service is 
equal to the distribution function of the busy period, for the given sys- 
tem. This is the curious and puzzling result mentioned in the abstract. 
The first result of this paper is the proof that, for Poisson input to a 
single server, the two distributions are alike only for exponential service. 
For more than one server, it is plausible that the same thing is true, 
though no easy line of proof seems open. For many servers, each with 
exponential service time, proof of identity may be given in a way parallel 
to the single server case. A second result is the determination of the delay 
distribution for Poisson input to a single server with constant service 
time; this result belongs in the book with that of Burke 3 on random 
service for the same system. 

II. LAST-COME FIRST-SERVED DELAY FOR POISSON INPUT TO A SINGLE 
SERVER 

For last-come first-served service, an arrival finding the server busy is 
delayed until all subsequent arrivals finding the server busy or just com- 
pleting service have been served. The waiting demands at the arrival 
epoch have no effect on this delay because the new arrival goes to the 
head of the line. 

Following Takacs, 4 the delay distribution may be derived as follows. 
Take the arrival rate as a, the service rate (reciprocal of the average 
service time) as b, the service time distribution function as B(t), and the 
delay distribution function as G*(t) —the star indicating that G*(t) 
differs from G(t), the distribution function for the busy period, in the 
starting epoch of the corresponding intervals. 

Note first that the (stationary) distribution function between an arbi- 
trary time epoch in a service interval and the epoch of next service 
completion is given byf 

C(() = b [ [1 - B(n)]du, (1) 



t The earliest proof of this result known to me is in Palm 8 ; it may have been 
proved much earlier. 
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so that 

y(s) = I e-'dCit) (2) 

= (6/a)[l - /3(a)], 



with 



0(a) = f c-'d.B(t). 
Jo 



Note that -y(a) = j8(a) if and only if /3(s) = 6(6 -f a)" 1 , or what is the 
same thing, £(/) = 1 - c - '''. 

Now consider the interval between the arrival epoch of a delayed 
demand and the epoch of the next service completion. The probability 
that this interval has length (y, y + dy) is dC(y). The probability of n 
arrivals in this interval, when its length is y, is the Poisson term 
e~ ay {ay) n /n\ If n = 0, the delay of the given arrival for last-come first- 
served service is simply ('(/); if n = 1, the delay consists of the interval 
to the first service completion phis the busy interval occasioned by this 
arrival and all subsequent arrivals during subsequent service periods, the 
distribution function for which is G(t). Hence for n = 1 the delay is the 
convolution of C(t) and G(t). In the same way, for n = 2 the delay is 
the convolution of ('( I ).(/( / 1 and (!{ t ), and so on. 

If G„(l) is the distribution function for the sum of n variables, each 
with distribution function G(t), and G (t) = 1, G^t) = G(t), then the 
conditional delay distribution function is given by 

G*(t) = f T, c-" ( W£ G n (t - y) dC(y) . (3) 

If 

T*(a) = / e~ 3l dG*(t) 

Ji\ 

and r(a) has a similar significance, then the transform of (3) ist 

T*(s) = y[s + a - aT(s)}, (4) 

with 7(a) defined by (2). Note that T (a) is determined by Takacs' equa- 
tion 

r(a) = P[s + a - aV(s)]. (5) 

t A similar and equivalent result appears in Wishart*; note that Wishart con- 
siders the unconditional delay distribution function, which is less directly com- 
parable (than the conditional) with the distribution function of the busy period. 
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Hence T*(s) = T(s) when and only when y(s) = /3(s). As already 
noted, this implies B(t) = 1 - e~~ ht , the exponential service distribution. 

The common distribution function for exponential service time has 
been given by Takacs, 4 and is obtained as follows. First, since /3(s) = 
b(b + s)~\ the solution of (5) is 

_, , _ a + b + s - V(q + b + s) 2 - 4ab 
USJ _ ~ 2* 



Next, the inverse of thisf is 
dx 



(7(0 = ( -^7=e- (a+b)x h(2xVab), P = a/b, 
Jo x\/p 



(6) 



with h{x) the Bessel function of the first kind and imaginary argument. 
An equivalent expression due to Vaulot, 1 which seems better adapted to 
numerical work, is 

1 _ G(t) = - f' tt e ' Abt si » 2 x > A = 1 + p - 2>/p cos x. (6a) 

III. DELAYS FOR POISSON INPUT TO MANY EXPONENTIAL SERVERS 

The modifications of the argument above for c fully accessible servers, 
each with exponential distribution of service time (exponential servers, 
for brevity), are relatively minor. First, in the derivation of the busy- 
period distribution, a service period in the single-server case is replaced 
by the interval between the epoch at which the last idle server becomes 
busy and the epoch at which the first of c busy servers becomes idle. 
The distribution function of this interval is that of the least of c random 
variables, each with the same exponential distribution, the service time 
distribution. If this function is H(t) and the service time distribution is 
B(t) = 1 - e' bt , then 

1 - H{t) = [1 - B(l)] e = e- bc '. (7) 

Hence the distribution function G c (t) is obtained from its single server 
(exponential service time) correspondent (7(0 simply by replacing b 

by be. 

The last-come first-served delay distribution is proved identical simply 
by the remark that H(t) is also the distribution of the interval between 
an arbitrary point in the interval for which H(t) is the distribution 
function and its termination, because H{t) is of exponential form. 

t Ref . 7, pair 55G.1, p. 59. 
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Fig. 1 shows u comparison of conditional delay curves for various orders 
of service at an occupancy p = a/bc of 0.9. The orders are order of ar- 
rival, random, and inverse order of arrival (last-come first-served). The 
abscissae are values of an auxiliary time variable u = bet. The ordinates 
arc values of F(u), the probability that the delay of a delayed arrival is 
at least u. The curve marked "envelope" is an upper bound for all pos- 
sible orders of service; actually it is a plot of 1 — v(u) with 

v(u) = (1 - p) J2p"v n (u) 

and v„(u) as denned in the introduction. 

A detailed study of delays for last-come first -served service and the 
busy period, with and without defections from waiting, is planned for a 
later paper. 
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Fig. 1 — Comparison of delay curves for various orders of service occupancy 0.0. 
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IV. DELAY FOR SINGLE SERVER WITH CONSTANT SERVICE TIME 

Although this is a limit case of (4), both delay and busy-period dis- 
tribution functions are obtained more easily, following Takacs for the 
latter, from the classification of the busy period by number served. 

Note first that the generating function for number of arrivals in an 
interval with distribution function (7(0 is given by 

= £ flV [" e- v (y n /nl) dC(y) (8) 

n=0 ''O 

= y(a — ax). 

The corresponding generating function for number of arrivals in a service 
interval is given by 

p(x) = j8(a - ax) (9) 

and, by (2), 

„.(*)_ I_ZJ!<£> P -a/b. (10) 

p(l - X) 



Note that 



P P 



pp,* = 1 - p — Pi — ' ' • — Pn ■ 

Now write /„* for the probability of n .services in a delay period for 
last-come first-served service, and 

./■*(.(•) = E/«v 

n = l 

for its generating function; /„ and f(x) are the corresponding entities for 
a busy period. Then, following Takacs, 

/■* = Pi7. , 

.A* = px*h + pffi, 

and 

/„* = Pl*/»-l + P2*E/if-l-i + P**Hfifkfn-l-i-* + •", (ID 

/*(*) = .rp*[/(,r)l = .r 7 [a - a/(.r)]. (12) 
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Note that the corresponding result for/(.r) obtained by Takacs is 

/(*) = xfla - af( ..-)]. (13) 

Hence, by (2), (12), and (13), 

rM-tn+mv > = a " K i,4> 

Note that 

f( x ) = f l*> - axfix)3'[a - af(x)]. 
x 

Nolo also that (14) may be rewritten as 

pf*(x) = 1 + (.r - 1)[1 -fix)]- 1 

= 1 + U - l)g(x), 

where 

g(x)[i -/(.*•)] = l, 

0(0) = g = 1. 
Hence 

pfn* = g.,-1 - g„ . (16) 

For constant service time with service rate b, 
B{t) = 0, t < b~\ 

= l, / > ir 1 

and I3(s) = e~' lb . Hence 

/(*) = xc-" + " f(r \ 
whose solution is given by 



/,.,:, = £ if^lll 



x" 



„=1 111 
so that 



(14a) 



(is; 



/ \n— 1 —no 
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Then, using (15), the coefficients g a of the generating function g(x), are 
given by 

go = l, 

0i = h = e" P , 

&2=/ 2 +/l 2 = (1 +PK 2P , 

g 3 = h + Wi +/i S = (1 + 2p + 3p 2 /2!)e- 3p . 
These results suggest writing 

From (15) it follows that 

?n = E/tf.-;- (18a) 

j-i 

From this and (17) and (18), it is found that 

^- = Z(t : )0-+l) J '"V H -y-u-i- (19) 



From (19) it is found in succession that g ll0 = 1, 0»i = w — 1, B 2 = 
n(w - 2), all of which are contained in the single formula g uk = n'~' 
(n — k). If this is correct, substitution in (19) shows that the following 
must be an identity: 

n k -\n - k) = t( k ) U ~ ^(n - ; - l)*-^(n - ft - 1). (20) 

.'=0 V/ 

Equation ( 20) is in fact one of the forms associated with Abel's generali- 
zation of the binomial formula; it is a special case of Equation (lb) of 
Salie, 8 which is due to Jensen 9 (of whom it is proper to remember here 
that he was chief engineer of the Copenhagen Telephone Company as 
well as a mathematician of note). 
Thus, finally, f 

g« = e- p £,(n ~ k)n k -'p k /k\ (21) 

k=0 

which, with (16), determines /„*( re). 

Note also that the polynomials g n = g„(p) appear also in the dual 

t An equivalent result has been found independently by Frank A. Haight, of 
the University of California. 
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problem of the busy period of a single server with exponential service 
time and regular arrivals considered by Connolly. 1 " 
The first service has the distribution function 

St(t) = U, t ^ b'\ 

(22) 
= 1, t £ 6" 1 

and the distribution function for n services is 

S n (t) = 0, f £ (n - 1 )/>-', 

= bt, (n - 1)6 _1 < t < nb~\ (23) 

™ 1, nb~ < t. 

So 

G*(t) = T,fn*S a (t) (24) 

n=l 

is the (conditional) delay distribution function for last-come first-served 
service with Poisson input to a single server with constant service time. 
The distribution function for the busy period (Takacs 4 ) is 

0(0 = f^"(p») H /"!, (25) 

with [bt\ the integral part of bt. 
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